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1 Introduction 



In these lectures I will show, through three examples, how current lattice calculations are 
able to tell us interesting things about the continuum physics of non-Abelian gauge theories. 

My first topic concerns the glueball spectrum. The physics question here is: where, in the 
experimentally determined hadron spectrum, are the glueballs hiding? I will first summarise 
what lattice calculations tell us about the continuum glueball spectrum of the SU(3) gauge 
theory. I will then discuss what this tells us about the masses of the corresponding 'bare' 
glueballs in QCD. I will then turn to the experimental spectrum with some discussion of 
the interpretation of the observed states in the quark model. Finally I will pinpoint the 
experimental states most likely to have large glueball components. 

The second topic concerns topological fluctuations in the 577(3) gauge theory. Here the 
simplest physics question is: are these fluctuations large enough to be consistent with the 
observed large 7/ mass? Thanks to Witten and Veneziano this is a question that can be 
posed in the pure gauge theory. We shall see that the fluctuations of the topological charge 
do indeed have the required magnitude. Along the way I will discuss the problems with 
topology on a lattice. I then move onto the much less straightforward question concerning 
the structure of these vacuum fluctuations: e.g. what is the size distribution of instantons? I 
will discuss some preliminary calculations that show the mean size to be about p ~ 0.5/m. I 
will also point to some intriguing evidence for a long distance polarisation of the topological 
fluctuations. 

My third topic concerns the physics of SU (N c ) gauge theories as a function of the number 
of colours, N c . This will be mainly in 2+1 dimensions, since that is where we have good 
calculations. I will describe calculations of the mass spectrum for N c < 4 which explicitly 
show that for N c >2 mass ratios are independent of N c up to a modest ~ V-^c correction. 
This is a very elegant result: it tells us that all the apparently different SU(N C ) theories are 
actually one single theory, SU (00), to a reasonable first approximation. There is some very 
preliminary evidence, as I will show, that the same is true for D — 3 + 1. 

2 Glueballs in QCD 

Ideally I should be telling you what happens when you simulate QCD with realistically light 
quarks. But it is going to be a few years yet before I can do that. What current lattice 
Monte Carlo calculations are able to provide is predictions for the low-lying mass spectrum 
of the continuum SU(3) gauge theory without quarks. These states are glueballs - there 
being nothing other than gluons in the theory. If you want hadrons with quarks then you 
can propagate quarks in this gluonic vacuum and then tie such propagators together so 
that the object propagating has the appropriate hadronic quantum numbers. That is to 
say, you calculate hadron masses in the relativistic valence quark approximation. (This is 
usually referred to as the 'quenched approximation' to QCD.) The spectrum one obtains this 
way is a remarkably good approximation to the observed hadron spectrum. This is not too 
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surprising: one reason we were able to learn of the existence of quarks in the first place is 
because the low-lying hadrons are in fact well described by a simple valence quark picture. 

Suppose that we begin with our SU(3) gauge theory and then couple to it 3 flavours of 
very heavy quarks. Initially the spectrum will contain the usual light glueball spectrum, 
supplemented by a spectrum of very heavy quarkonia that can be well accounted for in 
terms of a valence quark potential model. Let us now gradually reduce the quark masses 
towards their physical values. In principle the glueball and quarkonia states might entirely 
change their character once their masses become comparable. However, as we remarked 
above, the experimental light quark spectrum still seems to retain the essential features of 
valence quark physics. If the quarkonia are not qualitatively alterred, it seems reasonable 
to think that neither will the glueballs be. Of course if a quarkonium state and a glueball 
are close enough in mass they will mix. However there is reason to believe that this mixing 
is weak. The reason is the Zweig (OZI) rule: hadron decays where the initial quarks all 
have to annihilate are strongly suppressed. The classic example is the <fi meson. Such a 
decay may be thought of as quarks — > glue — ► quarks. Glueball mixing with quarks should 
therefore be y/OZl suppressed. As should glueball decays into hadrons composed of quarks. 
The existence of such a suppression is supported by a recent lattice calculation M . 

The picture we have in mind is therefore as follows. The glueballs will only be mildly affected 
by the presence of light quarks. They will decay into, say, pions but their decay width will be 
relatively small; and there will be a correspondingly small mass shift. Only if there happens 
to be a flavour singlet quarkonium state close by in mass will things be very different, because 
of the mixing of these nearly degenerate states. In this context we expect 'close by' to mean 
within ~ lOOMeV. So we view the glueballs in the pure SU(3) gauge theory as being the 
'bare' glueballs of QCD which may mix with nearby quarkonia to produce the hadrons that 
are observed in experiments. All this is an assumption of course, albeit a reasonable one. If 
true it tells us that the glueballs, whether mixed with quarkonia or not, should lie close to 
the masses they have in the gauge theory. So we now turn to the calculation of those masses. 

I shall begin by reviewing the available lattice calculations, placing a particular emphasis on 
exposing the sources of systematic error in arriving at a final mass prediction in MeV units. 
I do this in some detail, so that you are able to judge for yourselves the credibility of lattice 
mass estimates. 

As we shall soon see, it is only for three states that the lattice calculations are reliable 
enough that we can extract continuum predictions: the lightest scalar, ++ , tensor, 2 ++ , and 
pseudoscalar, 0~ + , glueballs. For other states we do have calculations for one or two values 
of the lattice spacing a, but that is not enough to extrapolate to a = 0. Nonetheless the 
lattice results strongly suggest that glueballs with other J PC are heavier f2|. 

The lightest glueball is the ++ , and it is for this state that we have the most accurate 
lattice predictions. Although there have been recent estimates for the mass that appear to 
differ, e.g. 1.55 ± 0.05GeV § and 1.74 ± 0.07GeV |, |, we shall see that this difference 
is illusory. In fact the apparent difference reflects different ways of extrapolating to the 
continuum limit and different ways of introducing physical MeV units into the pure gauge 
theory. That is to say, it reflects particular systematic errors which we need to estimate and 
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this we try to do. Since we find that the various lattice calculations are consistent, we are 
able to carry out a global analysis that provides the best available glueball mass estimate. 
We find m ++ = 1606 ± 73 ± 130MeV where the first error is statistical and the second is 
systematic. Performing a similar analysis for the lightest tensor and pseudoscalar glueballs 
we find m 2 ++ = 2200 ± 73 ± ISOMeV and m -+ = 2100 ± 73 ± 130MeV. 

From the experimental and phenomenological point of view the scalar sector is complex, and 
I will review the current state of play. As we shall see, there appear to be too many scalar 
states to be explained as quarkonia, and the strongest candidates for states with large gluonic 
components are the / (1500) and the /(j =0 ?)(1710). This possibility is strongly reinforced 
by the fact that these are the only experimental states that are compatible with the lattice 
mass estimate in the previous paragraph. 

In order to assess which states have the largest gluonic components, it is important to 
understand the mixing between nearby quarkonia and glueballs. I shall briefly discuss what 
happens in the case of mixing between a glueball and the lightest singlet and octet scalar 
quarkonium multiplets. 

For this review I have drawn very heavily on 01]. 
2.1 Calculating glueball masses 

I begin by briefly reminding you of some general aspects of lattice calculations. (For more 
detail see the lectures by Chris Michael at this School and the books by Creutz || and 
Montvay and Munster ||.) In a lattice calculation (Euclidean) space-time is discretised, 
usually onto a hypercubic lattice. What we calculate is the mass spectrum of the discretised 
theory, but what we actually want is the corresponding spectrum of the continuum theory. 
Because the theory is renormalisable, the effects of the lattice spacing, a, on physical length 
scales will vanish as a — > 0. Since QCD has effectively one length scale, ~ one 
expects the effects of the discretisation to become negligible once a « lfm. The same 
is true of the gauge theory without quarks. (By '1/m' I really mean some characteristic 
physical length scale. How one introduces actual fermi units into the pure gauge theory is 
something we shall return to below.) The lattice spacing is varied by changing the value of 
the bare (inverse) coupling which appears in the lattice action: (3 = Q/g 2 . Since the theory 
is asymptotically free, we know that in order to approach the continuum limit, a — > 0, we 
need to take g 2 — * and so (3 — > oo. Indeed for sufficiently small g 2 one can determine the 
relationship between a and g 2 in low-order perturbation theory. (Although in practice the 
latter does not work well for the range of couplings currently accessible.) Reducing a makes 
the calculation numerically more intensive for various reasons. An obvious one is that if we 
wish to maintain a constant volume, the number of lattice sites grows oc 1/a 4 . It is only 
in recent years that calculations for very small values of a have become practical. As we 
shall see, the mass calculations I shall use here have been performed over a range of lattice 
spacings 0.18/m > a > 0.06/m, corresponding to couplings 5.7 < (3 < 6.4. 

We now outline the main steps in a lattice Monte Carlo calculation of a glueball mass with 
a view to exposing the main different sources of systematic error. (So we ignore various 
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unilluminating technicalities and skate over various qualifications that are irrelevant for this 
purpose.) 

• Extracting lattice masses. The first step is to calculate the glueball spectrum on a lattice, 
for a given space-time volume V and for a particular lattice spacing a. The Euclidean time 
translation operator is e~ Ht , where H is the Hamiltonian of the theory. Thus the correlation 
function of an operator with some particular J PC quantum numbers will, for large enough 
values of t, vary as e~ mot where mo is the lightest mass with those quantum numbers. 
We calculate such propagators numerically for t large enough that we see this asymptotic 
exponential decay. From the exponent we extract the mass. Since t is given in lattice units 
(t = an where n is the number of lattice spacings) the exponent is m Q t = m an and so 
what we actually obtain is the mass in lattice units, i.e. am$. The statistical errors are 
straightforward to estimate. However there is also a systematic error that has to do with 
determining the range of t where the asymptotic exponential dominates. We do not attempt 
to quantify this error but simply note that it will become increasingly important for the 
heavier glueball states (such as the tensor and pseudoscalar) where the exponential decrease 
with t of the 'signal', and hence its immersion into the statistical noise, is more rapid. 

• Finite V corrections. The second step is to determine the corrections due to the fact 
that the volume is finite. For most quantities the functional forms of the leading \arge-V 
corrections are known theoretically; typically they will be of the form 5m oc exp{— camL} 
where am and L are the mass-gap and lattice size in lattice units |7]] . By doing calculations for 
a variety of volumes at some chosen value of a, any unknown constants in these expressions 
can be fitted and the resulting formulae can be used to apply corrections at other values of a. 
Since the calculations we use here have been performed on periodic volumes of typical sizes 
1.5 to 2.0 fm, these corrections are smaller than our typical statistical errors; but it is hard 
to determine them more accurately than that. This provides another source of systematic 
error. 

• Finite a corrections. The calculated glueball mass will depend on a; and it will do so in 
two ways. Firstly, it is obtained in lattice units, am. This trivial dependence is removed 
when we take the ratio of two masses: am\jam2 = m\jm2- The non-trivial a dependence is 
due to the distortion of the dynamics by the discretisation. For the lattice action we use and 
for quantities such as glueball mass ratios it is known that the leading small-a correction 
is 0(a 2 ) H . So for small enough a we can extrapolate our calculated mass ratios to a = 



where m may be chosen to be mi or m.2 or some other physical mass: the difference be- 
tween these choices is clearly higher order in a 2 . Such neglected higher order terms in the 
extrapolation are another source of systematic error. 

• Introducing the MeV scale. Having obtained the continuum mass spectrum in the form of 
mass ratios, we want to express the masses in usual MeV units. This can be done if at least 
one of the masses corresponds to a quantity whose value is known in MeV. For example 
the potential between heavy quarks is linear for large separations: V(r) ~ or where a is 
called the string tension. In simple string pictures for high J hadrons a is related to the 
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slope, a', of Regge trajectories by a' = l/2na and this provides the conventional estimate 
y/o ~ 400 — A50MeV ||. Knowing the continuum value of mc/y/a we can now express 
ma in MeV units. Of course one may distrust this particular argument for the value of a. 
An alternative is to calculate quark propagators in the pure gauge theory and from these to 
form hadron propagators. From the asymptotic exponential decay of the latter we obtain 
quarkonium masses in the (relativistic) valence quark approximation. For example, we can 
obtain the continuum limit of the mass ratio m p j\/a. Setting m p = 770MeV we obtain a 
value for y/a and hence for the glueball mass mg. Or we might set the scale using the <fi- 
meson or the nucleon instead of the p. While the mass spectrum one obtains in the quenched 
approximation is remarkably close to that which is experimentally observed, it is not exactly 
the same and so these different ways of setting the physical scale will lead to slightly different 
glueball masses. This is a source of systematic error. 

The first three types of systematic error can be made arbitrarily small by sufficiently im- 
proving the numerical calculations (in obvious ways). The error in setting the MeV scale is 
qualitatively different. It is intrinsic to working within the quenched approximation. How 
do we estimate it? There are some quantities which we know are going to be sensitive to the 
absence of vacuum qq fluctuations, e.g. the rf mass or the topological susceptibility. These 
should obviously not be used to set the MeV scale in the pure gauge theory. There are other 
quantities which we expect to be no more sensitive to the absence of vacuum qq fluctuations 
than the glueball masses themselves and which are therefore suitable quantities with which 
to attempt to set the MeV scale. These include the masses of typical quarkonia such as the 
p meson, the <fi meson, the nucleon, certain matrix elements, etc. Of course the ratios of 
these quantities cannot be exactly the same in full and quenched QCD and so the scale we 
extract will vary according to which of these quantities we choose to use. The extent of this 
variation can be used as an estimate of the systematic error. 



2.2 The lightest glueballs 

The values of the glueball masses that we shall use are from [JTT], [TD|, |^, and those of the 



string tension are taken from |T(], 0, |T^|. As we have mentionned already, only the ++ , 



2 ++ and 0~ + glueballs are determined accurately enough that a continuum extrapolation is 
possible. In J| you can find mass estimates for glueballs of widely varying J PC , obtained for 
a very small lattice spacing. These do suggest that the three masses we shall obtain are in 
fact the lightest ones. There is clearly an urgent need for a new generation of lattice glueball 
calculations which will provide information on a much larger part of the continuum mass 
spectrum. For a preview of what these are likely to look like, see ]T3| . 



The first step is to take ratios of masses so that the scale, a, in which they are expressed 
cancels. We choose to take ratios of the glueball masses, amc, to a^/a since this latter 
quantity has been very accurately calculated. Now we remarked above that for small enough 
a the leading discretisation effects in such mass ratios are 0(a 2 ). So for small enough a we 
expect 

m G (a) m G (a = 0) 2 

~~r~T\ = ~T( nT + ca a ( 2 ) 
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We take all the available mass values and try to fit them using eqn |2|. If a good fit is not 
possible we assume that this is because the largest values of a used is too large for the 0(a 2 ) 
correction to be adequate. So we drop the mass corresponding to the largest value of a and 
try again. We keep doing this until we get a good fit. We find that the ++ glueball can be 
well fitted in this way over a range of lattice spacings 0.18/m > a > 0.06/m, corresponding 
to couplings 5.7 < (3 < 6.4. (Note that where we employ fermi units, these have been 
introduced using y/o~ = A40MeV; a value that will be made plausible later on.) Such lattice 
spacings are small enough that we are not surprised that higher order terms in a 2 should be 
small. We show the mass ratios in Fig [l] together with the best fit of the form in eqn We 
obtain the continuum mass ratio: 

^ = 3.65 ±0.11 (3) 

The fit is obviously a very good one (with a confidence level of 85%) and it is clear that the 
calculations of the different groups are entirely consistent with each other. 

In Fig |2] and Fig |3] I show corresponding plots for the 2 ++ and 0~ + glueballs. The former is 
well determined, and we obtain the continuum ratio 

^ = 5.15 ±0.21. (4) 

It is however clear that our control over the _+ is marginal. This translates into a very 
large error when we perform the continuum extrapolation: 

^ = 4.97 ±0.58. (5) 

We now wish to transform the above continuum glueball masses to physical MeV units. 
There are several reasonable ways to do this and how they differ will give us an estimate of 
the systematic error intrinsic to introducing physical units into a theory that is not quite 
physical. 

As we remarked earlier, one way is to infer from the observed Regge slopes that ^Ja ~ 
400 — A50MeV. This estimate does however suffer from being somewhat model dependent. 
An alternative is to take lattice calculations of the mass of the p and extrapolate the ratio 
m p/V a t° the continuum limit, just as we did for the glueball. The only difference with 
eqn |2| is that the leading correction will be O(a) rather than 0(a 2 ). We do this for two recent 
'state-of-the-art' calculations: the GF11 collaboration JliJ and the UKQCD collaboration 
fl~5| who use quite different discretisations for the quark action. UKQCD uses an improved 
action which should have smaller discretisation errors. We plot the two sets of ratios in Fig ^ 
with their corresponding continuum extrapolations. These give us 



and 



p - 1.72 ±0.08 GF11 (6) 



—f- = 1.78 ± 0.08 UKQCD (7) 
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respectively. It is reassuring that these two calculations are entirely consistent in the con- 
tinuum limit, despite the fact that they have very different lattice discretisation corrections. 
(Indeed one can argue that the dominant correction in the UKQCD calculation will be 0(a 2 ) 
rather than 0(a).) To extract a value for a in MeV units we average the above results, set 
m p = 770 MeV" and so obtain sjo = 440±15MeV. We observe that this is entirely consistent 
with the scale we inferred from Regge slopes; but the argument is much cleaner here. 

The error on a is largely statistical. We now need to estimate the systematic errors as well. 
These are discussed in detail in ||. There we consider errors due to finite volume corrections, 
to extrapolations in the quark mass, to uncertainties in our estimates of the lattice values 
of a and to using the K* or nucleon to set the scale rather than the p. We estimate a ±8% 
systematic error in total. This leads to our final estimate for the value of the string tension 
as being: 

yja = 440 ± 15 ± 35MeV (8) 
where the first error is statistical and the second is systematic. 

We can now use this value in eqns 0- 1^ to express our glueball masses in MeV units. We 
obtain 

m 0++ = 1.61 ±0.07 ±0.13GeV (9) 
m 2++ = 2.26 ±0.12 ±0.18GeV (10) 

and 

m -+ = 2.19 ±0.26 ±0.18GeV. (11) 

The first error combines the statistical errors on the glueball and p continuum extrapolations, 
and the second is our estimate of the systematic error. This, then, is our best lattice 
prediction for the lightest glueballs prior to any mixing with nearby quarkonium states. 

Where are the quarkonia? We are interested in flavour-singlet states because those are the 
ones that can mix with glueballs. It would clearly be very useful if lattice calculations were to 
provide estimates for the masses of the uu + dd and ss mesons prior to their mixing with glue. 
One could then introduce the mixing with the scalar glueball using the (standard) formalism 
described below and compare the resulting states with the experimental spectrum. There are 



some indications from recent lattice calculations |T6|, [L7| that there is a scalar ss state close 
to the scalar glueball mass. However these quenched calculations do not try to incorporate 
the essential quark annihilation contributions and so must be regarded as indicative at best. 
We will therefore turn now to experiment and phenomenology. 



2.3 Experiment and phenomenology 

An illuminating (if arbitrary) starting point is provided by first considering the lightest 2 ++ 
mesons. In the quark model [Tj| these are very similar to the scalars: the quark spins are 
aligned in both cases and the only difference is that the net spin is parallel to the unit 
orbital angular momentum for the tensors and antiparallel for the scalars - which, in quark 
potential models, leads to minor differences in the masses Hl9| . However in the real world we 



expect the tensors to have narrower decay widths than the scalars because their decays into 
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light pseudoscalars require non-zero angular momentum and corresponding near-threshold 
suppression factors. Thus they should be easier to identify experimentally - and that is their 
interest for us here. 



Experimentally we find the following lightest tensor states f20|| . There is an isoscalar / 2 (1270) 
with width, T ~ 185MeV, a second isoscalar /2(1525) (r ~ 76MeV), an isovector a2(1320) 
(r ~ 107MeV) and a strange isodoublet K*(U30) (T ~ WOMeV). The first isoscalar 
decays mainly into pions while the second decays mainly into strange mesons. Thus it is 
natural to infer that the /2(1270) is mainly uu + dd while the /2(1525) is mainly ss. We 
have a clear nonet of tensor mesons and we note that the splittings are exactly what one 
would expect from a mass- difference between strange and non-strange (constituent) quarks 
of ~ lOOMeV. Thus the lightest tensors provide a nice illustration of the quark model at 
its most successful: the mesons fall into 577(3) multiplets with a modest symmetry breaking 
driven by the m s — m n mass difference. (For convenience we shall adopt the shorthand 
notation nn for -^{uu + dd) in the following.) 

We now turn to the ++ mesons. The easiest such meson to see experimentally should 
be the strange isodoublet and the lightest such state turns out to be the i\~g(1430) with 
T ~ 290MeV. Note that this is close in mass to the corresponding tensor although, as 
expected, its decay width is much larger. (This is what makes the scalars very much harder 
to identify experimentally.) There is also a candidate isovector a (1450) (r ~ 270MeV). 
In addition there are two isoscalars, the /o(1370) (r ~ 300 — 500MeV) and the /o(1500) 
(r ~ 12QMeV). The former decays mainly into pions and so one would suppose it to be 
composed mainly of non-strange quarks. 

So far all this looks much like the tensor nonet. However there is a puzzle. The /q(1500) 
does not have the obvious decays for a predominantly ss state. Moreover if we compare 
its width to the other scalars then we see that it is remarkably narrow - particularly for 
a state with, apparently, a large non-strange component. This motivates us to look for 
other nearby scalar states. There is evidence that the /j(1710) is, or contains, a ++ state. 
Moreover the predominant 2-body decays of this state involve strange quarks (KK or rjrj). 
The state is relatively narrow T ~ 175MeV. So if we include this state into our discussion, 
we no longer have an obvious ss problem. However now we have three isoscalar states in 
the 1370 — 1710MeV mass region - too many for a quark model nonet! Since, as we have 
seen, lattice calculations predict a scalar glueball in the 1600 MeV mass region, with a width 
comparable to that which we observe for the /o(1500) and /j(1710), it is natural to conjecture 
that the three observed isoscalars are in fact the results of mixing between two quarkonium 
isoscalars and the lightest scalar glueball. This is the scenario explored in [^3], |2l], |22| with 
some differing assumptions. 

All this might be a convincing picture if it were not for the existence of some lighter scalar 
states that we have so far failed to mention. These are the isoscalar / (980) and the isovector 
a (980). Both are narrow (r ~ 40 — lOOMeV). (There may also be an extremely wide iso- 
iscalar in the 400 — 1200MeV mass range.) Since there is no nearby strange isodoublet, 
these states certainly do not fit into the usual quark picture where mesons should fall into 
approximate SU(3) multiplets with modest symmetry breakings driven by the m s — m n 



mass difference. In fact these mysterious states were interpreted some time ago as 
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being loosely-bound KK molecules - recall that they are narrow and occur very close to the 
KK threshold. If we accept this interpretation - and it has been widely accepted as being 
plausible - then we can ignore these states for our purposes and the interpretation of the 
previous paragraph remains. I should stress that alternative interpretations of this spectrum 
do exist, but I will ignore these here and instead refer you to |4| for a detailed discussion. 

The picture we are thus led to is one where the / (1370), the / (1500) and the / (1710) 
are the result of mixing between the glueball and the would-be ss and nn scalar quarkonia. 
We argued above that these quarkonia belong to the same nonet so the ss will be heavier 
than the nn. Moreover the relative strengths of the glueball mixings will be determined on 
symmetry grounds. 

What are the constraints on the mixing? 

First the output masses of the mixing should correspond to the / (1370), / (1500) and 
/o(1710). The first state is very broad and so we shall allow its mass to lie in the region 
Mi e [1.31, 1.40](j?eV with a preference for 1.37GeV. The second state is both quite narrow 
and well-defined and so we fix its mass to M 2 = l.bGeV. The third state appears to be 
again quite narrow but the precise location of the ++ component is still quite uncertain. 
We shall consider the range M 3 e [1-64, 1.80]GeV with a preference for l.HGeV . 

There are also some constraints on the input parameters. Given the lattice predictions, 
a generous range for the glueball mass would be mc £ [1-40, 1.80]GeV. We also have a 
qualitative constraint on the glueball-quarkonium matrix element of the Hamiltonian: it 
should be small, ~ 0(100) MeV or less. Similarly we expect a glueball to be narrow, and 
so any state with a large glueball component should be narrower than one would otherwise 
expect. As for the quarkonia, we expect m sS > m nn , with a mass difference in the ~ 150 — 
200MeV ball-park. And the mixing matrix elements should be close to their SU '(3) flavour 
values. There are further, and important, constraints that arise from the observed decays of 
the three output states and again I refer to ||] for a more detailed discussion. 

Here I shall give one example of what appears to be an acceptable mixing scheme. We assume 
that prior to mixing what we have are the 'bare' quarkonia with masses m sS = l.QIGeV and 
m nn = 1.36GeV and a bare glueball with mass mg = 1.48GeV. These are not eigenstates 
of H and so the mass matrix will not be diagonal. We can write the latter as 

/ m G z 

z m sS 
V y/2z 

Winn I 

(The factor of a/2 follows from the fact that ^{< uu\+ < dd\}H\G >= y/2 < uu\H\G >.) 
The physical masses, after mixing, are the eigenvalues of this mass matrix. With an accept- 
ably weak mixing of z — 62MeV we obtain as output masses 1.31, 1.50, 1.64GeV for the 
/o(1370), /o(1500) and / (1710), which is perfectly acceptable. The overlaps of the output 
states can be written in terms if the bare input states as follows: 



i 


fin 


fis 


fiG 


/o(1370) 


0.87 


0.10 


-0.49 


/o(1500) 


0.48 


-0.43 


0.77 


/o(1710) 


0.13 


0.90 


0.42 
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We see that most of the glueball resides in the /o(1500) and the /q(1710) is mainly an ss 
quarkonium. Such a mixing scheme is qualitatively compatible with the observed decays of 
the output states. This and other possible mixing scenarios are discussed in M For some 



other recent discussions of mixing in this context see [21, 22 



2.4 Conclusion 

We now know quite accurately the lightest scalar glueball mass in the pure gauge theory 
in units of the string tension. If we translate this into physical units, as described in this 
lecture, we find that the glueball should appear around ~ 1.60 ± 0.15GeV. This is just 
the mass range where one naively expects the two scalar flavour singlet quarkonia to lie. If 
they do then all these states will inevitably mix to some extent, even though we expect the 
dynamical mixing parameter to be weak. It is therefore intriguing that experimentally not 
only are there definitely states in this mass range, the / (1370) and the / (1500), but there 
is evidence for a third, the /o(1710). We presented a sample mixing scheme in which the 



glueball mainly resides in the /o(1500). This is not a new suggestion | 21f . While it is still 
a little too early to come to a convincing conclusion, it certainly seems that it will not be 
very long before we are able to do so. 

I have emphasised the scalar glueball because that is where most of the recent interest 
has been. What about the tensor and pseudoscalar? In looking for glueball candidates 
experimentalists naturally look for states that appear in 'gluon-rich' processes. For example 
in J 1 1\) decays, or in Pomeron-Pomeron collisions, or states that have large rjr] decays. This 
picks out the /o(1500) and /o(1710) amongst the scalars. Amongst the tensors this picks 
out the /2(1900) and the G(2150). These are, of course, in the right ball-park for the 2 ++ 
glueball as predicted by lattice calculations: ~ 2.26 ± 0.22GeV. So here too, things look 
interesting. 

With the pseudoscalar, on the other hand, we have little reason to feel smug. First, the 
lattice calculations are poor: a mass estimate of ~ 2. 16 ±0.32 barely qualifies as an estimate 
at all. Moreover the obvious experimental candidate is the t(1490) seen in J/ip decays. We 
should however remember that this state is special: it has the quantum numbers of the 
vacuum topological charge. Which brings me smoothly to my next topic. 



3 Topological fluctuations 

As you know, SU(N) gauge fields in 3+1 dimensions possess a topological charge. PS[ . And 
it is the topological fluctuations of the gauge fields that are the reason why the rf has a 



mass ~ lGeV rather than being almost a Goldstone boson |26|]. Moreover there is good 
reason to think that these fluctuations lead to the spontaneous breaking of chiral symmetry. 
The reason is that isolated instantons produce zero modes in the Dirac operator; these mix 
with each other, and shift away from zero, when the instantons are not isolated (as in the 
real vacuum). It is not hard to imagine that this might leave a non-zero density of modes 
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close to zero, and this would suffice to break chiral symmetry (via the Banks-Casher formula 



1 27]). This occurs explicitly in some instanton models (see Shuryak's lectures at this School 



and [28. and has been observed in some lattice calculations p5 |. It is also possible that 
instantons may affect the properties of some hadrons, as the near-zero modes may lead to 
large contributions to the valence-quark propagators. This is more speculative. All this to 
say that topology is interesting. In addition it is intrinsically non-perturbative. Here I will 
tell you some things that we have learned by studying topology in lattice gauge theory. 

The first thing I need to address is the fact that when we discretise space-time we lose 
topology in a formal sense, since any field on a discrete set of points can be smoothly 
deformed to a trivial field. One should not get too excited about this; the same happens 
with dimensional regularisation. If the space-time dimension is not exactly 4 we have no 
topological winding. Since the theory is renormalisable we expect that, as the cut-off is 
removed, a — > 0, we recover all the properties of the continuum theory. We shall see how 
this occurs for topology in the lattice gauge theory. 

So in this lecture I will address the following topics. First I discuss the basic ambiguity with 
defining topology on a lattice; and show why it does not really matter. Then I discuss the 
practical problem of calculating the topological charge of fields that have fluctuations on all 
length scales. This is a large subject and I will focus simply on the one technique that goes 
by the name of 'cooling'. I will then move onto the first bit of physics: the Witten-Veneziano 
fomula that relates the strength of the topological fluctuations in the pure gauge theory to 
the mass of the rf in QCD. Finally I will attempt to give some insight into the structure of 
the topological fluctuations in the vacuum. 



3.1 A basic ambiguity and why it does not matter 

Let me start by giving an explicit example. A lattice gauge field is defined by a set of SU(2) 
matrices on the links, I, of the lattice: {U{\. (I will stick here to SU(2) for simplicity.) 
Consider now the following continuum gauge potential for an instanton of size p centered at 
x = 0: 2 

K(x) = -^— 2 g-\x)d,g{x) (12) 
x z + p z 



with 



Xq + IXjCTj . . 

9( x ) = ~ a/2 ( 13 ) 



Let us now translate this field by a/2 in each direction so that it is centered at x 
(a/2, a/2, a/2, a/2), i.e at the centre of a lattice hypercube. Define a lattice field by: 



x+afi 



Ufa)=Vexp J Al(x)dx (14) 

X 

For p ^> a the lattice discretisation is very fine compared to the instanton core in which the 
instanton action and topological charge density, Q(x) = j^e^ vp aTr{F^ v (x)F pa (x)} , reside. 
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In this case any reasonable definition of the topological charge (see below for an example) 
will assign a topological charge of Q = 1 to this lattice field. Suppose we now continuously 
reduce p. The lattice field will also vary continuously Eventually we will have p <C a. At 
this stage the instanton core, which is at the centre of a lattice hypercube, will be very far, 
in units of its size, from any of the lattice links. Thus even the nearest link matrices on the 
lattice will be arbitrarily close to pure gauge. That is to say, the field configuration will be 
the same as that due to a gauge singularity located at x — x. In this case any reasonable 
definition of the topological charge will assign a topological charge of Q = to this lattice 
field. Thus we have passed continuously from a field with Q — 1 to one with Q = 0. 

Does this raise a fundamental problem with topology on the lattice? The answer is: not 
really. At least if we are interested, as here, with the limit a — > 0. Let me argue why this is 
so. 

Consider the density of topological fluctuations as a function of their size p. This is only an 
unambiguous notion if p <C ^ where £<j is the typical dynamical length scale of the theory 
(e.g. 1 fermi in QCD or the inverse mass gap in the pure gauge theory). In that case we 
know the density of these 'instantons': 

do 1 8tt 2 

D(p)dp = — — e (15) 
p p 4 

where the represent factors varying weakly with p. You recognise in this equation the 
scale-invariant integration measure; also a factor to account for the fact that a ball of volume 
p 4 can be placed in 1/p 4 different ways in a unit volume; and finally a factor arising from the 
classical instanton action, Sj = 8n 2 /g 2 , with perturbative fluctuations promoting the bare 
g 2 to a running g 2 (p) in the usual way. This last sentence is the most important one for us 
here: if we substitute for g 2 (p) in eqn |l^ we find that 

D(p)dp oc p 6 dp : SU(3) (16) 

with a power p 7//3 in the case of SU(2). So , because of the scale anomaly, the number of 
instantons rapidly vanishes as p — ► rather than diverging as 1/p 5 . 

On the lattice this density will change as follows if a -C p and a -C 

D(p)^D(p)x{l + 0(^)} (17) 

P 2 

Now, suppose a lattice field configuration is to be smoothly deformed from Q = 1 to Q = 0. 
This requires a topological fluctuation to be squuezed out of the lattice, as described above. 
While we do not know much about the structure of the original fluctuation (it will typically 
be on a size scale ~ ^ which is beyond the reach of our analytic techniques) we do know that 
if the lattice spacing is sufficiently small then to reach p ~ a the 'instanton' will have to pass 
through sizes £^ ^> p ^> a. In this region the density is calculable as we saw above, with a 
probability that is very strongly suppressed; at least as ~ (p/£d) 6 for SU(3). So the changing 
of Q is conditional upon the involvement of field configurations whose probability— > as 
a — > 0. Thus, as we approach the continuum limit this lattice ambiguity vanishes very 
rapidly. We knew that this had to happen because the theory is renormalisable and the 
lattice is surely a good regulator. It is, however, nice to see it happen explicitly. 
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3.2 Cooling 



Having convinced you that it makes sense to discuss topology on the lattice, I now want to 
discuss how you can calculate it. The obvious ways are three. 

• Calculate the zero- modes of the Dirac operator /D[A|. Q will equal the difference between 
the number of left and right handed zero-modes. 

• Interpolate a smooth gauge potential. If Q ^ and if space-time is compact (as it is here: 
a hypertorus) then we will need more than one patch. The net winding of the transition 
functions between the patches equals Q. 

• We calculate the topological charge density Q{x) = -^je llU p rj Tr{F IJL y{x)Fp (J {x)}. Then 
Q = JQ(x)d 4 x. 

All these methods have been explored in lattice calculations. I will focus on the last because 
it is the simplest and because it immediately tells us something about the size and location 
of the core of the 'instanton' (since that is where Q(x) is localised). 

We need a lattice operator that becomes Q(x) in the continuum limit. This is easy. Recall 
(see Chris Michael's lectures) that if we define the plaquette matrix, Up U (x), as the ordered 
product of link matrices around the corresponding plaquette of the lattice, then it is easy to 
see that U^ u (x) = 1 + a 2 F^ u (x) + .... and hence that ||T 



Ql(x) = ^e^ pa Tr{Up V {x)U pa {x)} = a 4 Q(x) + 0(a 6 ). (18) 

If we apply this formula to an instanton of size p (discretised as described above) then we 
find, as expected, that Ql = J Qi{x)dx = 1 + 0(a 2 /p 2 ). 

Applied to the real vacuum however Ql{x) has problems. The operator is dimensionless 
and so 0(a 6 ) actually means terms like ~ a 6 ~ a 6 FD 2 F etc. For smooth fields these 
are indeed 0(a 6 ). However realistic fields (those that contribute to the path integral) have 
fluctuations all the way up to frequencies of 0(l/a). So if we take matrix elements of Ql(x) 
in the vacuum the contribution of the high frequency modes to the 0(a 6 ) terms will be 
8Ql(x) ~ a 6 x 1/a 6 ~ 0(a°). In practice this contribution is suppressed by some powers of ft 
that can be calculated in perturbation theory (since these contributions are short- distance). 
Thus in the real world Ql{x) possesses interesting topological contributions that are of 
order a 4 oc e~ c/3 (since g 2 (a) oc l/log(aA)) and uninteresting ultraviolet contributions that 
are oc l/ft n . So as we approach the continuum limit, ft — > oo, the latter dominate and we 
are in trouble. 

Actually things are a little worse than this. Like other composite lattice operators, Ql{x) 
possesses a multiplicative lattice renormalisation factor: Zq ~ 1 + d j ft + ... P^| . This looks 
innocuous, and indeed in the continuum limit it obviously is. However it turns out that the 
value of d is such that Zq <^ 1 precisely in the range of values of ft where current lattice 
calculations are performed. 

There are different ways to deal with these problems. I will describe a particularly simple 
technique ||33|| . The idea rests on the observation that the problems are all caused by the 
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ultraviolet fluctuations on wavelengths ~ a. By contrast, if we are close to the continuum 
limit, the topology is on wavelengths /)> a. One can therefore imagine taking the lattice 
fields and locally smoothing them over distances ^> a but <C p. Such a smoothing would 
erase the unwanted ultraviolet fluctuations while not significantly disturbing the physical 
topological charge fluctuations. One could then apply the operator Ql{x) to these 'cooled' 
fields to reveal the topological charge distribution of the vacuum. 

How do we cool a lattice gauge field? The simplest procedure is to take the field and generate 
from it a new field by the standard Monte Carlo heat bath algorithm subject to one crucial 
modification: we always choose the new link matrix to locally minimise the plaquette action. 
Since TrU^x) measures the variations of the link matrices over a distance a, minimising 
the plaquette action is a very efficient way to erase the ultraviolet fluctuations. (Obviously 
there are many possible variations on this theme.) 

Thus the notion is that we take our ensemble of N gauge fields, {U I=1, " ,N }, perform a 
suitable number of cooling sweeps on each one of these, so obtaining a corresponding ensemble 
{Uc =1 '"' N } of cooled fields, and then extract the desired topological properties from these 
cooled fields. As one might expect there are ambiguities for realistic values of a. As we cool, 
topological charges of opposite sign will gradually annihilate. This changes the topological 
charge density but not the total value of Q. Eventually this leads to a very dilute gas of 
instantons. As we cool even further these isolated instantons will gradually shrink and will 
eventually shrink within a hypercube and at this point even Q will change. (This is for a 
plaquette action on a large enough volume: other actions may have other effects.) Of course 
when an instanton becomes narrow it has a very peaked charge density and is impossible 
to miss. So we certainly know when it disappears out of the lattice and can, if we think 
it appropriate, correct for that. All this to say that cooling is a good way to calculate the 
total topological charge, but is not necessarily a reliable way to learn about the topological 
charge density. 

Let me give you an example of how it works. I have produced a sequence of 90 Monte Carlo 
sweeps on a 12 4 lattice at j3 = 2.5. These field configurations are separated by just one 
heat bath sweep (i.e. each link matrix has been changed only once) so one expects the long 
distance physics on neighbouring field configurations to be almost identical. That is, the 
value of Q should change little. In Fig. || I show you the values of Ql when calculated on 
these fields. They jump all over the place and are nowhere near the expected integer values. 
Now let us cool each of these 90 configurations with 25 cooling sweeps. Calculating on 
these one finds a dramatic difference - as shown in Fig. [5[ (The cooled charges differ from 
integers because there is an 0(a 2 /p 2 ) error which is substantial for our not-so-small choice of 



We now have a technique for calculating the topological charge of a lattice field. As we 
have seen, its application requires some care. Another reason to take care is the following. 
The probability of instantons with p ~ a depends on the lattice discretisation. If the lattice 
action is much less than the continuum action then even when a — > there may be a finite 
density (per unit physical volume) of these lattice artifacts. They might survive a couple 
of cooling sweeps. For realistic values of a the gap in scales between a and ^ is not very 
large and so such a narrow instanton, if sitting on a background field due to a physical but 
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moderately narrow anti-instanton, might expand under cooling rather than shrinking. So 
it might occasionally survive our 20, or whatever, cooling sweeps and bias our calculations. 
Obviously this latter effect disappears as a — > 0. Another correction is due to the fact 
that with finite a one loses the tail of the continuum density that extends to p < a. This 
correction also disappears as a — > 0. To deal with these problems one needs to perform a 
suitable scaling analysis. For example suppose we define a topological charge Ql{.Pc) which 
only includes topological charges for which p > p c = an c . As a — > Ql{Pc) should become 
independent of p c for a growing range n™ m < n c < n™ ax : with an™ in large enough to exclude 
any p ~ a artifacts and n™ ax growing exponentially with (3. We shall demonstrate a simple 
calculation of this kind next. 



3.3 The topological susceptibility 

Since < Q >= (we have no ^-term) the simplest quantity we can calculate is < Q 2 >. 



However there is a much better reason for calculating it: it is directly related |34|, |35[ to the 
rj mass: 

,2 o2 



« S ^^2^ ( < + <-2<) (19) 
If we put in the experimental numbers into the above formula, then we get 

Xt ~ (180 MeVf (20) 

In this relation the susceptibility, Xt, is that of the pure gauge theory - which is something 
that we can calculate. We now describe such a calculation and see what happens. 

The values of < Q 2 > that I am going to use come from calculations that I have done in 
the past; for SU(3) they come from [J7| while for SU(2) they come from 



All these calculations have been performed using 20 to 25 cooling sweeps. The charge Ql is 
then calculated on these cooled lattice fields. This charge is non-integer because the smallest 
instanton charges suffer significant 0(a 2 /p 2 ) corrections. However such small instantons have 
very peaked densities and are easy to identify. One can then estimate the corrections due to 
these small instantons, and shift Ql to the appropriate integer topological charge Q. This 
has been done in all these calculations. 

In Fig. ^| I plot Xt/ \J& against the string tension, a 2 a in lattice units. We expect the leading 
lattice corrections to this dimensionless mass ratio to be 0(a 2 /p 2 ), as in eqn |l|. So we would 
attempt a continuum extrapolation of the form 

Xt(a) Xt(0) , „ 2 



y/<r(a) y/a(0) 



+ ca 2 a (21) 



which is a simple straight line on the plot. As we see the calculated values are consistent 
with this functional form. In addition to the susceptibility calculated from the total charge 
we also calculate Q with narrow instantons removed. The particular cut we have used is to 
remove all charges whose peak density is greater than l/16n 2 . This corresponds to removing 
instantons with p < 3a. In doing this we are largely removing lattice artifacts with p ~ a 
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which, because of their environment, survive the cooling. Of course we are also removing 
a part of the small-p tail of the continuum density, D(p). So we are not saying that this 
is necessarily a better measure of the continuum susceptibility. What we do wish to check 
is that both these susceptibilities are consistent when extrapolated to the continuum limit. 
And as we see in Fig. |] this is indeed so. We extract: 

lim ^- = 0.437 ±0.020 ±0.015 (22) 

Here the first error is statistical and the second is a systematic error estimated using the 
two different extrapolations. If we now plug in our favourite value for the string tension, 
yjo = 440MeV, we obtain 

Xt = (192 ± YIMeVf (23) 

Of course we should really incorporate the uncertainty in the value of y/a, most of that 
coming from setting MeV units, and this would increase the error to ±20MeV. Irrespective 
of such details it is clear that the pure gauge theory susceptibility is indeed consistent with 
being large enough to drive the large rf mass. 

In Fig. [7| I show the corresponding plot for the SU(2) susceptibility. Again we see that the 
two susceptibilities we calculate are consistent when extrapolated to the continuum limit. 
Note that the difference between the two at corresponding values of a is much greater for 
SU (2) than for SU (3). This is because the SU (N) running coupling, g 2 (p), runs much faster 
to zero as N j and so any finite-a ambiguities between physical and ultraviolet topological 
charges rapidly decrease. We extract: 

lim ^- = 0.470 ±0.011 ±0.015 : 577(2) (24) 

which translates to 

Xt = (207 ± 9MeV) 4 : SU{2) (25) 

if we use \Ja — 440MeV. Unfortunately we do not know the r/ etc. masses for the SU(2) 
theory so we cannot say if this is what is expected. And neither does it make much sense 
to introduce MeV units by using the 577(3) value of sjo in this way. It does however allow 



me to compare to two recent estimates of 200 ± 15MeV jyj and of 230 ± 30MeV |42| which 
have been obtained using quite different methods, but a similar MeV scale. As we see, the 
values are all consistent - which is reassuring. 



3.4 Vacuum topological structure 

The calculation of Q on the lattice is relatively straightforward. If we continuously deform 
a continuum gauge field so as to minimise the action then we necessarily get driven to the 
semiclassical multi-instanton minimum. We do not change topological charge sectors since 
these are separated by infinite action barriers. Cooling is just a naive lattice version of this 
procedure and so it should work in this way for a — > 0. As we have seen, it does indeed seem 
to work well. 
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However there are other things we would like to know about the instanton density if we are 
interested in gauging the possible influence of instantons on chiral symmetry breaking, the 
quarkonium spectrum and related physics. For example a picture in which instantons have 
a mean size of ~ 1/3 fm and are relatively dilute produces a plethora of interesting physics 
p8| . Can we say something about such more detailed features of the vacuum topological 
structure? 

This question involves both technical and conceptual ambiguities. There are two main tech- 
nical problems. The first is that as we cool a gauge field its topological structure changes: 
instantons will change their sizes and nearby instantons and anti-instantons will annihilate. 
These quantities are not invariant under minimising the action; they are not even quasi- 
stable. So we must do as little cooling as possible and trust only those features that we 
find to be relatively insensitive to cooling. If the vacuum that one obtains is a dilute gas 
of instantons then there is no further difficulty However what one might expect, and in- 
deed finds, is a relatively dense gas of overlapping charges. This raises a difficult problem 
of pattern recognition, particularly for the larger instantons whose density, Ql{x) will be 
very small (it obviously varies as 1/p 4 ) and which are most likely to overlap with several 
other large instantons. Such large instantons will, in any case, not be completely smooth, 
and may possess multiple peaks ('ripples') since we are trying to minimise the number of 
cooling sweeps used. These problems have been addressed in different ways by several recent 
calculations of this kind HTI, f4"2L IP], I think it is fair to say that all these should 



be regarded as exploratory. However, to whet your appetite I am going to show you some 
results from an analysis that I have been involved in ]4"5|] . What I am not going to do here 
is to tell you anything at all about the pattern recognition algorithms etc. 

I mentionned a conceptual difficulty. Instantons are semiclassical objects. The real vacuum 
has fluctuations on all scales and there is no reason to think that the topological charge 
resides in instanton-like objects. Does it even make sense to talk of a distribution in p? 
This raises all kinds of interesting questions, including: do we care? After all, models such 
as that in [28j would not claim to be more than simplifications that are appropriate to the 



physics being considered. Perhaps we should regard a minimal amount of smoothening as 
performing such a simplification upon the fluctuating gauge fields? I am not going to address 
these questions any further here, but you should be aware of their existence. 



Let me now move to some results of this kind of calculation IE5I. The calculations are in the 



pure SU(3) gauge theory. We have used lattice fields that have been generated and stored 
(for other purposes) by the UKQCD Collaboration. We are performing calculations on 16 3 48 
and 32 3 64 lattices at /3 = 6.0, 24 3 48 lattices at (3 = 6.2, and 32 3 64 lattices at (3 = 6.4. The 
point of the 2 lattice sizes at (3 = 6.0 is to check for finite volume effects - especially for the 
very large instantons. The various lattice spacings enable us to check whether the features 
we identify have the right scaling properties to survive into the continuum limit. In addition 
we do the calculations for various different numbers of cooling sweeps so that we can check 
whether these features are insensitive to cooling or not. 

In Fig. |] I show the size distribution, -D(p), with the size p expressed in units of l/y/o- = 
1/y/K. The peak is around p ~ 0.5 fm. This is not very sensitive to (3 or to the lattice size 
or to the number of cooling sweeps. The total number of charges is of course sensitive to the 
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amount of cooling. We do, however, note that unless we go to a very large number of cools, 
the vacuum is dense with a great deal of overlap between the charges. 

It is not straightforward to test for scaling, since a cooling sweep is intrinsically non-scaling. 
What we do find is that if we tune the number of cooling sweeps so that the total number 
of charges is independent of j3, then the detailed densities, D(p), seem to scale as well. 

Two other aspects of D(p) are of particular interest: the fall-offs at small and large p. The 
former is interesting because it provides a check on our calculations: we know that D(p) oc p 6 
for a <C p <C 1/\A>"- Of course we can only approach these severe inequalities by looking for 
the trend as we increase (3. We do indeed find a tendency to approach something close to 
this functional behaviour. (Which, in any case, is modified by powers of log p.) By contrast 
large p is interesting because we do not know what happens to large instantons in a confining 
vacuum. If we used the semiclassical formula with a coupling that froze at large p then we 
would get -D(p) oc 1/p 5 as a — ► oo. What we find is a much stronger suppression: 

D(p) oc : p > ±- (26) 



Before closing this subject let me move to a quite different and more subtle aspect of the 
vacuum structure. In a dilute gas the correlation between the sign of a charge and the sign 
of the total charge will be independent of size p. In Fig. |9] I show the correlation we actually 
find. There is a striking size-dependence. Charges smaller than average tend to have the 
same sign as Q, larger charges tend to have the opposite charge. Since the sign of Q simply 
tells you which sign wins out on that configuration, this tells us that net charge of the small 
instantons wins out over the net charge of the larger instantons. This suggests a picture 
where the very large instantons are polarised, so that the small instantons sitting on this 
background tend to have the opposite charge thoughout the volume. Since the charge of the 
small instantons wins out, this means that the large instantons are actually overscreened. 
Of course this would be a lot to glean from the one plot; however this intriguing picture is in 
fact supported by more detailed calculations. It leads to quite striking effects. Suppose for 
example we decide to calculate the total charge that includes instantons smaller than some 
value p c : Q p < Pc - The fluctuations are shown in Fig. [H]. We note that if we were to limit 
ourselves to instantons with p < l/^/cr then we would obtain a topological susceptibility 
that is ~ 10 times greater than the total susceptibility! It is clearly important to investigate 
the effect of these structures on quark propagators and hence on the observable physics. 



3.5 Conclusions 

In this lecture I hope that I have convinced you that simulations of lattice gauge theories 
can be both useful and interesting in telling us something about the topological fluctuations 
of the vacuum. Some of what I have shown you, particularly concerning the long distance 
polarisation of the vacuum, must be regarded as preliminary - to a potentially embarrassing 
degree. The calculations of the topological susceptibility, on the other hand, are now quite 
reliable. They show us that the topological vacuum fluctuations are indeed large enough to 
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drive the large rf mass. The theoretical expectations with which we compared our result 
arise, most straightforwardly, from arguments about what happens in SU(N) QCD at large 
N. This takes me smoothly to the final of my three topics in these lectures. 



4 SU(N C ) gauge theories for all N c 

Quantum Chromodynamics is an 577(3) gauge theory coupled to 3 lightish quark colour 
triplets. If we change the gauge group to 577(2), 577(4), 577(5), ... then we obtain an infinite 
set of, a priori, quite different theories. However from an analysis of Fenyman diagrams to all 
orders ||B one finds that the N c — > oo limit of such SU(N C ) theories is smooth if we vary the 
coupling as g 2 oc 1/N C . This suggests that it should be possible to describe SU(N C ) gauge 
theories as perturbations in powers of l/N c around SU (oo) PEfl , at least for large enough N c . 



Moreover if one assumes confinement for all N c , then one can easily show [fllj, [|7| that the 
phenomenology of the SU (oo) quark-gluon theory is strikingly similar to that of (the non- 
baryonic sector of) QCD. This makes it conceivable that the physically interesting SU(3) 
theory could be largely understood by solving the much simpler 5£/(oo) theory. If all the 
SU (N c ) theories down to SU (3) can be treated in this way, then this represents an elegant 
and enormous theoretical simplification. 



I do not have the time to review []48fl this subject here, but let me at least indicate something 
of what is involved, albeit using arguments that lack any rigour. The constraint that g 2 oc 
1/N C is easy to motivate. Consider inserting a gluon loop into a gluon propagator. We have 
added two triple-gluon vertices and this gives a factor g 2 . At the same time the sum over 
colour in the loop gives a factor of N c . (Not N 2 because the colour of the incoming/outgoing 
gluon is fixed.) So we have a total factor ~ g 2 N c . Now such loops can be inserted any 
number of times, so if we want a smooth large- N c limit ( at least in all-order perturbation 
theory) then we clearly need to impose g 2 oc 1/N C . Assume therefore that we have done so. 
Consider a typical meson decay, e.g. p — > mr. This requires the production of a qq pair. So 
the decay width contains a factor of g 2 oc 1/N C . But if we have confinement all the hadrons 
are colour singlets and so we do not acquire any compensating factors from summing over 
the colours of the decay products. Thus the decay width is suppressed by a factor 1/N C . 
That is, at large N c hadrons do not decay. This is like a 'narrow-width' caricature of the 
real world. In fact this is a reasonable first approximation to the hadrons we know; mostly 
their widths are very much smaller than their masses. A similar argument tells us that 
mixings, e.g. between mesons and glueballs, are suppressed. This is reminiscent of the OZI 
rule discussed in my first lecture. This (and much more) suggests that the SU (oo) theory is 
indeed a first approximation to SU(3). 

All this motivates us to try and answer by explicit calculation some basic questions: 

• does a non-perturbative calculation support the (all-orders) perturbative argument for a 
smooth N c —>■ oo limit? 

• is N c —>■ oo confining? 
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• does such a limit really require g 2 oc 1/N C \ and if so what does this mean when we have a 
running coupling? 

• is SU(N C ) ~ SU (oo) only for iV c 1 or is it the case down to N c = 3, or even N c = 2? 

• what is the SU(oo) mass spectrum? 

There have been a number of interesting computational explorations of the lattice SU(oo) 



theory (for a review see based on the fact that it can be re-expressed as a single plaquette 
theory p0[ . Unfortunately this scheme makes no statement about the size of the leading 
corrections to the N c = oo limit, and so gives us no clue as to how close N c = 3 is to N c = oo. 

In this lecture I will describe an extremely straightforward approach to this problem. I will 
simply calculate the properties of SU(N C ) gauge theories for several values of N c and so 
determine explicitly how the physics varies as N c increases. I will only look at the pure 
gauge theory, but there are good reasons for believing that the inclusion of quarks will not 
alter any of our conclusions (except in some obvious ways). Ideally I would like to present 
you with accurate calculations in 4 dimensions. Unfortunately, at present what I can provide 
you there is very rough and tentative. But I will make up for that by describing the results 
of the corresponding, but much more precise, calculation in 3 dimensions. Of course you will 
want to know what is the relevance of such a 'substitution'. I will come to that shortly. 

My study was originally motivated by the observation that the C = + sector of the light 
mass spectrum turned out to be quite similar in the D = 2 + 1 SU(2) |Tl], |52|] and SU (3) 



53fl theories. (This also appears to be the case in D = 3 + 1, although there the comparison 
is weakened by the much larger errors.) One reason for this might be that both are close to 
the N c — oo limit. In that case we would have an economical understanding of the spectra 
of SU (N c ) gauge theories for all N c : there is a common spectrum with small corrections. 

A second, more practical reason for studying N c —>■ oo was my interest in obtaining some 
model understanding of the structure of glueballs. Models are of interest even if they are 
very approximate in comparison with the results of simulations. A good model will embody 
the essential degrees of freedom in a problem and show how this leads to the main features 
of the physics, within some transparent and plausible approximation scheme. A model could 
provide the intuition necessary for an economical understanding of the role of glueballs in 
a wide range of contexts. For example, if one understands the structure of glueballs, one 
can make crude but reliable estimates of glueball-quarkonium mixing, of glueball decay, of 
the effects of dynamical quarks etc. In many cases the approximations made in the model 
include the neglect of decays and mixing. These are features of the SU(oo) rather than of 
the SU(3) theory, and so it would be better to test the model against the spectrum of the 
former theory. One example that has been of particular interest to me is the flux tube model 
of glueballs [53], The formulation of this model is identical for all iV c > 2. However, 
because the model does not incorporate the effects of glueball decay, it should presumably be 
tested against the iV c — > oo spectrum since it is only in that limit that there are no decays. 
It is also the case that many theoretical approaches are simpler in that limit. An example is 
provided by the recent progress in calculating the large N c mass spectrum using light-front 
quantisation techniques [ pwjj . 
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The D = 2 + 1 analysis that I shall present here is based on my calculations over the last 
few years of the properties of SU (N c ) gauge theories with N c = 2,3,4 and 5. In D = 3 + 1 
what I have done is to perform some SU(4) calculations to supplement what is known about 
577(2) and S77(3). My strategy is the very simple one of directly calculating the mass 
spectra of these theories and seeing whether they are approximately independent of N c . The 
calculations are performed through the Monte Carlo simulation of the corresponding lattice 
theories, using the standard plaquette action. In the D — 2 + 1 case the calculations are 
very accurate and we are able to extrapolate our mass ratios to the continuum limit prior to 
the comparison. In the 3+1 dimensional case our 577(4) calculations are not good enough 
for that, and our comparisons with SU(2) and 577(3) are correspondingly less precise. Some 
of the SU{2) results have been published [pj], [5^] as have brief summaries of the results 



discussed here [pR [59|]. A long paper is in preparation. 



4.1 D=2+l ~ D=3+l ? 

While one might naively expect that the D — 2 + 1 and D — 3 + 1 gauge theories would be 
so different as to make a unified treatment misleading, this is not in fact so. Theoretically 
the D = 2 + 1 theory shares with its D = 3 + 1 homologue four important properties. 

• Both theories become free at short distances. In 3 dimensions the coupling, g 2 , has di- 
mensions of mass so that the effective dimensionless expansion parameter on a scale I will 
be 

gl{l) =lg 2 ^0 (27) 
In 4 dimensions the coupling is dimensionless and runs in a way we are all familar with: 

In both cases the interactions vanish as / — > 0, although they do so much faster in the 
super-renormalisable D = 2 + 1 case than in the merely asymptotically free D = 3 + 1 case. 

• Both theories become strongly coupled at large distances. This we see immediately by let- 
ting I | in the above formulae. Thus in both cases the interesting physics is nonperturbative. 

• In both theories the coupling sets the mass scale. In 3 dimensions it does so explicitly: 

rrii = c l9 2 (29) 

In 4 dimensions it does so through the phenomenon of dimensional transmutation: the 
classical scale invariance is anomalous, the coupling runs and this introduces a mass scale 
through the rate at which it runs: 

rrii = qA (30) 

• Both theories confine with a linear potential. This is not something that we can prove 
by a simple argument. However lattice simulations provide convincing evidence that this is 
indeed the case. Note that although the D = 2 + 1 Coulomb potential is already confining, 
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this is a weak logarithmic confinement, Vc(r) ~ g 2 ln(r), which has nothing to do with the 
nonperturbative linear potential, V(r) ~ err, that one finds at large r. 

In addition to these theoretical similarities, the calculated spectra also show some striking 
similarities. 

• In both theories the lightest glueball is the scalar ++ with a similar mass m ++ ~ 4^/(7. 
In the C = + sector, the 2 ++ is the next lightest glueball (ignoring any excited scalars) with 
m 2 ++/m ++ ~ 3/2 in both cases. 

All this motivates us to believe that a unified treatment makes sense. 
Of course there are significant differences as well. For example: 

o There are no instantons in D = 2 + 1 non-Abelian gauge theories. This probably implies 
quite different quark physics. 

o The rotation group is Abelian. 

o The details of the mass spectrum are very different in 3 and 4 dimensions. A particularly 
striking difference is that in 3 dimensions the spectrum exhibits parity doubling for states 
with non-zero angular momentum: 

m J+ = mj- if J 7^ (31) 

This follows from the fact that angular momentum flips sign under parity - which can be 
defined in D = 3 as (x, y) — > (—x, y). The proof is elementary and I leave it as an exercise. 
Actually to show this you need the continuum rotation group. If you only have 7r/2 rotations 
then you find that you lose parity doubling for the J = 2 state. This can happen either 
because you are too close to the strong coupling limit or because the volume is too small and 
the rotation symmetry is broken by the boundary conditions. In either case this is a useful 
check: that we are effectively in the continuum limit and in an effectively infinite volume. 



4.2 2+1 dimensions. 

The calculations in 3 dimensions are performed in the same way as in 4 dimensions except 
that the computational problem is much more manageable. (Lattices grow as L 3 rather than 
as L 4 .) The basic steps are just as outlined in my first lecture, and you will have to trust 
me that all the checks have been performed sufficiently carefully! After extrapolating to 
the continuum limit I obtain the string tension and a mass spectrum in units of g 2 . Only 
the lightest portion of the mass spectrum is calculated, but that includes J PC states for 
J = 0, 1, 2 and P = ± and C = ± as well as one or two further excited states in many cases. 
I do this for SU(2), SU(3), SU(A) and SU(5). In the case of SU(2) there is no C = - sector. 
In fact the main point of the SU(5) calculation was to have 3 values of N c for the C = — 
states, so as to provide some control over their N c dependence. The SU(5) calculation is 
recent and has not been incorporated into the analysis that follows. 

I begin with the string tension, a, since it turns out to be our most accurately calculated 
physical quantity. We use smeared Polyakov loops J3SJ, to obtain a 2 a for several values of 



23 



the lattice spacing a. We then extrapolate the lattice results, using the asymptotic relation 
(3 = 4:/ag 2 , to obtain the continuum string tension in units of g 2 : 



y^L = i im JL a i a (32) 

The results for 577(2), 577(3) and 577(4) in D = 2 + 1 are shown in Table [I] and are 
plotted in Fig. [II]. We immediately see that there is an approximate linear rise with N c and 
we find that we can obtain a good fit with 

0.1974(12,^-™ (33) 

We obtain a similar behaviour with the light glueball masses (see below). 
Some observations. 

• For large N c , eqn |33| tells us that y/cr oc g 2 N c . That is to say, the overall mass scale of the 
theory, call it /i, is proportional to g 2 N c . In other words, in units of the mass scale of the 
theory 

s 2 K k < 34 > 

While this coincides with the usual expectation based on an analysis of Feynman diagrams, 
we note that here the argument is fully non-perturbative. 

• The string tension is non-zero for all N c and, in particular, for N c — > oo (when expressed in 
units of g 2 N c or the lightest glueball masses - see below). This confirms the basic assumption 
that needs to be made in 4 dimensions in order to extract the usual phenomenology of the 
large- N c theory. 



• In the pure gauge sector one expects (again from an analysis of Feynman diagrams) 
that the first correction to the large- N c limit will be 0(1 /N 2 ) relative to the leading term. 
The fit in eqn |33] is indeed of this form. We note that if we try a fit with a 0(1/N C ) correction 
instead (which would be appropriate if we had quarks) then we obtain an unacceptably poor 
X 2 (corresponding to a confidence level of only ~ 2% in contrast to the ~ 45% we obtain for 
the quadratic correction). We may regard this as providing some non-perturbative support 
for this diagram-based expectation. 

• The coefficient of the correction term in eqn ^ is comparable to that of the leading term, 
suggesting an expansion in powers of 1/N C that is rapidly convergent. Indeed one has to go 
to N c = 1 before the correction term becomes comparable to the leading term. While the 
577(1) theory is completely trivial, we note that the £7(1) theory has a zero string tension 
(in the sense that yjcr/g 2 = in the continuum limit). 

Let me now turn from the string tension to the mass spectrum. Recall that since we are in 
D = 2 + 1, states of opposite parity are degenerate as long as J ^ 0. This degeneracy is 
broken by lattice spacing and finite volume corrections. I will present the results separately 
for the P = + and P = — states so as to provide an explicit check on the presence of any 
such unwanted corrections. 
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I begin with the C = + spectrum since the SU(2) spectrum does not contain C = — states. 
In this case we have masses for three values of N c , and so can check how good is a fit of 



the kind in eqn 33. In Fig. 12. I plot the ratio mc/g N c against 1/N£ for a selection of 



the lightest states, G. On this plot a fit of the form in eqn |33]. will be a straight line and I 
show the best such fits. As we can see, the data is consistent with such a 1/N 2 correction 
being dominant for N c > 2. However what is really striking is the lack of any apparent iV c 
dependence for the lightest ++ and 2 ++ states. 

In Table [| I present the results of fitting the C = + states to the form 

m G R ^ + (35) 



g 2 N c °° iV c 2 

where i?^ = -pf*r- j . (Note that the errors on the slope and intercept are highly corre- 
lated.) The confidence levels of the fits are quite acceptable suggesting once again that for 
N c > 2 a moderately sized correction of the form ~ I/TV 2 is all that is needed. Note that 
since the variation with iV c is small, the exact form of the correction used will not have a 
large impact on the extrapolation to N c = oo (except in estimating the errors). 

These calculations confirm my earlier claim that the physical mass scale at large N c is g 2 N c . 
So if we consider ratios of m G to y/a (as was explicitly done in p9| ) we will find that they 



have finite non-zero limits as N c — > oo : that is to say, the large- N c theory possesses linear 
confinement. 

For the C = — states we only have masses for 2 values of N c and we cannot therefore check 
whether a fit of the form in eqn [35| is statistically favoured or not. However given that such 
a fit has proved accurate for the C = + masses and for the string tension down to N c = 2 
it seems entirely reasonable to assume that it will be appropriate for iV c > 3 for the C = — 
masses. Assuming this we obtain the results shown in Table [] for the N c — oo limit and 
for the coefficient of the first correction. The 'lever arm' on this extrapolation is, of course, 
shorter than for the C = + states and that leads to correspondingly larger errors. This will 
be dramatically improved once I incorporate the 577(5) spectrum into the analysis. 

The results in the Tables provide us not only with values for the various mass ratios in the 
limit N c — > oo but also, when inserted into eqn [33|, predictions for all values of 7V C . 

Finally, I should remark that I have also calculated the deconfining temperature, T c , for 
SU{2) 0| and for SU(3) |j53| . Extrapolating as in eqn |33], we find 

T c . . 0.079(23) . , 

0.1745(52) + (36) 



g 2 N c ' iV c 2 

Of course, extrapolating from N c = 2, 3 is much less reliable than extrapolating from N c 
3, 4, and so this relation should be treated with some caution. 



4.3 3+1 dimensions. 

Our knowledge of 4 dimensional gauge theories is much less precise. As far as continuum 
properties are concerned, quantities that are known with reasonable accuracy include the 
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string tension, the lightest scalar and tensor glueballs, the deconfining temperature and the 
topological susceptibility. As in 3 dimensions, the SU (2) and SU (3) values are within ~ 20% 
of each other, which encourages us to investigate the 577(4) theory so as to see whether we 
are indeed 'close' to N c = oo. Of course these £77(4) calculations are much slower than in 



D = 2 + 1 and the results I will present here are of a very preliminary nature [5c 



I use the standard plaquette action, and so the first potential hurdle is the presence of the 
well-known bulk transition that occurs as we increase (3 from strong towards weak coupling. 
To locate this transition I have performed a scan on a 10 4 lattice and found that it occurred 
at (3 = 10.4 ± 0.1. This corresponds to a rather large value of the lattice spacing, a, and so 
does not lie in the range of couplings within which we shall be working, i.e. (3 =10.7,10.9 
and 11.1. 

Our calculation consists of 4000,6000 and 3000 sweeps on 10 4 ,12 4 and 16 4 lattices at (3 =10.7, 
10.9 and 11.1 respectively. Every fifth sweep we calculated correlations of (smeared) gluonic 
loops and from these we extracted the string tension and the masses of the lightest ++ 
and 2 ++ particles, using standard techniques [|38| , 1Q| . These are presented in Table |j. We 



also calculated the topological susceptibility, a A Xt- The charge Q was obtained using the 
cooling method discussed in the previous lecture. These calculations were performed every 
50 sweeps. Overall this corresponds to rather small statistics and the errors are therefore 
unlikely to be very reliable. 

We see from Table (| that the most accurate physical quantity in our calculations is the string 
tension, a. Can we learn from it how g 2 varies with N c , just as we did in D = 2 + 1? We 
focus on a particular embodiment of this question: if we compare different SU(N C ) theories 
at a value of a which is the same in physical units, i.e. for which a^/a is the same, does the 
bare coupling vary as 1/N C , i.e. does (3 = 2N c /g 2 oc N 2 7 We perform this comparison for 
/?4 = 10.9, 11.1. (For convenience we shall label (3 by the value of N c , i.e. we write it as (3n c -) 
To find the corresponding values of j3 in SU(2) and SU(3) we simply interpolate between the 
values provided in (for example) [^, |l(J . Doing so we find that the values of (3 corresponding 
to /3 4 = 10.9,11.1 are (3 3 ~ 5.972(18), 6.071(24) and (3 2 ~ 2.442(9), 2.485(11) respectively. 
If we simply scale /3 4 by N 2 then what we would expect to obtain is (3 3 ~ 6.131,6.244 and 
02 — 2.725,2.775 respectively. Superficially the numbers look to be in the right ballpark, 
but in fact the agreement is poor. For example /3 2 = 2.725 and (3 2 = 2.442 correspond to 
values of a\Jo that differ by about a factor of 3. 

This disagreement should not, however, be taken too seriously, since it is well-known that 
the lattice bare coupling is a very poor perturbative expansion parameter. It is known that 
one can get a much better expansion parameter if one uses instead the mean-field improved 
coupling, gf, obtained from g 2 by dividing it by the average plaquette, < -^-TrU p > [|BT 



Defining (3* Nc = 2N c /gj(a) we find that f3 4 = 10.9,11.1 correspond to f3[ = 6.215,6.474 
respectively. Scaling (3{ by N 2 we would expect the equivalent S77(3) and SU (2) couplings 
to be given by /3| = 3.496,3.642 and = 1.554,1.619. What we actually find is that 
the equivalent couplings are /3| ~ 3.527(22), 3.649(28) and f3 T 2 = 1.561(10), 1.613(12). The 
agreement is now excellent. That is to say, if the SU(N C ) mean-field improved bare-coupling 
is defined on a length scale that is related to the physical length scale (y/cr) by some constant 
factor, then it varies as g 2 oc 1/N C . This is, of course, the usual diagram-based expectation. 
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In Fig. [13] I plot the scalar and tensor glueball masses, in units of y/a, as a function of N c . 
For N c = 2,3 we have used the continuum values. For N c = 4 the calculations are not 
precise enough to permit an extrapolation to the continuum limit and so we simply present 
the values that we obtained at (3 = 10.9 and 11.1. (We do not use the (3 = 10.7 values since 
they have large errors and there is the danger that the scalar mass may be reduced by its 
proximity to the critical point at the end of the bulk transition line.) Although the N c = 4 
errors are quite large, it certainly seems that there is little variation with N c for N c >2 and 
any dependence appears to be consistent with being given by a simple 1/N% correction. The 
fact that these mass ratios appear to have finite non-zero limits, implies that the large- N c 
theory is confining. 



we 



As mentioned earlier we have also calculated the topological susceptibility. In Fig. [14 
plot the dimensionless ratio Xt^/V 17 as a function of N c . Once again the N c =2 and 3 
values are continuum extrapolations of lattice values, while in the case of SU(4) we simply 
display the lattice values obtained at (3=10.9 and 11.1. As we remarked earlier one expects, 
semiclassically, very few small instantons for SU(4) and this is confirmed in our cooling 
calculations. This has the advantage that the lattice ambiguities that arise when instantons 
are not much larger than a are reduced as compared to SU(3), and dramatically reduced as 
compared to SU{2). This implies that the interesting large- N c physics of topology (and the 
related meson physics) should be straightforward to study. 



4.4 Conclusions. 

I have shown you the mass spectra and string tensions of gauge theories with N c = 2,3,4 
in 3 dimensions. We saw that there is only a small variation with N c and that this can 
be accurately described by a modest 0(1/N%) correction. That is to say, such theories are 
close to their N c = oo limit for all values of iV c > 2. We also saw that the large- N c theory 
is confining and that g 2 oc 1/N C when expressed in physical units. This confirms, in a 
fully non-perturbative way, expectations arrived at from analyses of Feynman diagrams. It 
simultaneously provides a unified understanding of all our SU(N C ) theories in terms of just 
the one theory, SU(oo), with modest corrections to it. In practical terms this means that, 
from the parameters in our Tables, we know the corresponding masses for all values of N c . 

Our calculation in 4 dimensions, while quite preliminary, does suggest that the situation may 
be much the same there. This is obviously something that needs to be done much better. 
Of particular interest are topological fluctuations at large N c : the instanton size density, 
D(p), loses its small p tail and lattice topology should become unambiguous. At the same 
time it would be interesting to investigate the behaviour of hadrons composed of quarks: for 
example the rj' mass should vanish as N c — > oo. There are many problems of real theoretical 
interest here. I hope some of you will choose to get involved. 
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N c 




2 


0.3350 (15) 


3 


0.5530 (20) 


4 


0.7564 (45) 



Table 1: The D = 2 + 1 SU(N C ) confining string tension. 



p 


Roo 


Rslope 




0.1974 ( 


12) 


-0.12 


( 1) 


0++ 


0.805 ( 


13) 


-0.06 


( 8) 


0++* 


1.245 ( 


27) 


-0.41 


(14) 


o-+ 


1.788 ( 


88) 


-0.48 


(56) 


2++ 


1.333 ( 


29) 


-0.08 


(18) 


2-+ 


1.340 ( 


40) 


-0.01 


(24) 


1++ 


1.946 ( 


75) 


-0.59 


(47) 


1-+ 


1.919 ( 


115) 


-0.18 


(75) 



Table 2: States with C = + in D 
the 1/N% correction in eqn 



2 + 1 : i?oo = lim 



m G 



iV c -*oo g 2 N c 



and R s iope is the coefficient of 



G 


Roo 


R s 


lope 


0~ 


1.18 ( 6) 


0.1 


(0.6) 


* 


1.47 (10) 


0.3 


(1.1) 


o+- 


1.98 (28) 


-0.4 


(2.7) 


2~ 


1.52 (14) 


0.9 


(1.4) 


2+- 


1.58 (13) 


-0.4 


(1.3) 


1— 


1.85 (15) 


-0.3 


(1.5) 


1+- 


1.78 (23) 


1.3 


(2.3) 



Table 3: As in Table 2 but for states with C — — . 



p 


a^Jo 


am ++ 


am 2 ++ 


10.7 
10.9 
11.1 


0.296 (14) 
0.229 ( 7) 
0.196 ( 7) 


0.98 (17) 
0.77 ( 8) 
0.78 ( 6) 


1.78 (34) 
1.20 (10) 
1.08 (10) 



Table 4: SU(4) in 4 dimensions; masses calculated at the values of (3 shown. 
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Figure 1: The scalar glueball mass: the GF11 values (x) and the rest 
extrapolation to the continuum limit is shown. 
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Figure 2: The tensor glueball mass: the GF11 values (x) and the rest 
extrapolation to the continuum limit is shown. 
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Figure 3: The pseudoscalar glueball mass. The best linear extrapolation to the continuum 
limit is shown. 
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Figure 4: The p mass: GF11 (x) and UKQCD (•) values. The best linear extrapolations to 
the respective continuum limits are shown. 
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Figure 5: The lattice topological charge: before cooling (x) and after cooling (•) the fields. 
Calculated on a sequence of fields separated by one Monte Carlo sweep. 
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Figure 6: The SU(3) topological susceptibility: from all peaks (•) and from peaks with 
Q{x P eak) < 1/167T 2 (x). Sample extrapolations to a common continuum limit are shown. 
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Figure 7: The SU(2) topological susceptibility: from all peaks (•) and from peaks with 
Q{x P eak) < 1/167T 2 (x). Sample extrapolations to a common continuum limit are shown. 
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Figure 8: The number of topological charges versus size at (5 — 6.2: after 23 (x), 32 (o) and 
46 (o) (under-relaxed) cooling sweeps. The size is in units of l/\/K where K is the string 
tension. 
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Figure 9: The average correlation between the sign of a charge of size p and the sign of the 
total charge of the gauge field; at f3 — 6.2 after 23 cooling sweeps 
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Q 2 (p<p cllt ) vs p cllt : (3 = 6.2: 23 sweeps 
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Figure 10: The average of Q 2 when only charges with sizes p < p c are included : at f3 — 6.2 
after 23 cooling sweeps 
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Figure 11: Continuum string tension versus number of colours in D 
in eqn KB. 



2 + 1. Line is the fit 
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Figure 12: Some continuum glueball masses, in D = 3, for 2,3,4 colours: ++ («), ++ *(x), 
2 ++ (*), _+ (o) and linear fits. 
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Figure 13: Lightest scalar (•) and tensor (o) glueball masses in D = 4. Continuum values 
for iV c = 2, 3 and lattice values ((3 = 10.9 and (3 = 11.1) for iV c = 4. 
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Figure 14: The topological susceptibility: continuum values for N c = 2, 3 and lattice values 
(P = 10.9 and = 11.1) for A^ c = 4. 
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